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Abstract 

We present a through discussion of motivations for and phenomenological issues 
in super symmetric models with minimal matter content and non-holomorphic soft- 
breaking terms. Using the unification of the gauge couplings and assuming SUSY is 
broken with non-standard soft terms, we provide semi- analytic solutions of the RGEs 
for low and high choices of tan(3 which can be used to study the phenomenology in 
detail. 

We also present a generic form of RGIs in mSUGRA framework which can be used 
to derive new relations in addition to those existing in the literature. Our results are 
mostly presented with respect to the conventional minimal supersymmetric model for 
ease of comparison. 
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1 Introduction 



Supersymmetry is an elegant symmetry for stabilizing the electroweak scale against strong 
ultraviolet sensitivity of the Higgs sector induced by quantum fluctuations. This symmetry, 
given that no experiment has yet observed any of the superpartners, cannot be operative 
at energies below the Fermi scale. This very constraint is saturated by breaking global 
supersymmetry explicitly via mass parameters 0{TeV) in such a way that the quadratic 
divergence of the Higgs sector is not regenerated. In more explicit terms, the action density 
of the minimal supersymmetric model (MSSM) which is based on the superpotential 

W ^ht tRQiHu + hb hnQjld + K fRLjld + IJ^HJId (1) 

as obtained after discarding all Yukawa couplings except those of the heaviest fermions, is 
augmented by additional terms (see, for instance, [1] for a review) 



htAttRQi^H^ + hMRQLHd + KArTRLLHd + /I'BH^Hd + ^ ^A„A„ + h.c. 



(2) 



which contain massive scalars, gauginos as well as a set of triscalar couplings among sfermions 
and Higgs bosons. The operators in (2) break supersymmetry in such a way that Higgs scalar 
sector does not develop any quadratic sensitivity to the UV scale. 

The soft-breaking terms in (2) do not necessarily represent the most general set of oper- 
ators. Indeed, one may consider, for instance, triscalar couplings with 'wrong' Higgs as well 
as bare Higgsino mass terms. Indeed, such terms have recently been shown to occur among 
flux-induced soft terms within intersecting brane models [2]. Historically, such terms have 
been classified as hard since they have the potential of regenerating the quadratic divergences 
[3]. However, this danger occurs only in theories with pure singlets, and in theories like the 
MSSM they are perfectly soft. Hence, the most general soft-breaking sector must include 
the operators 

l^'HuHd + htAtdRQLHl + HMrQlHI + KA'^trLlHI + h.c. (3) 

in addition to those in (2). Clearly, none of these operators mimics those contained in the 
superpotential (1): they are non-holomorphic soft- breaking operators. Note the structure 
of the triscalar couplings here; the triscalar couplings in (2) are modified by including the 
opposite- hypercharge Higgs doublet. 

In principle, the theory can contain both and ^' couplings. However, in what follows 
we will follow the viewpoint that the ii parameter is completely soft, that is, in the 
superpotential vanishes. This indeed can happen if the theory is invariant under global chiral 
symmetries [4] at high scale [5]. What is crucial about vanishing ^ is that it automatically 
solves the problem; the theory does not contain a supersymmetric mass parameter with 
a completely unknown scale. Indeed, in the MSSM stabilization of the parameter to the 
electroweak scale requires the introduction of gauge [6]- or non-gauge [7] extensions in which 
the vacuum expectation value (VEV) of an MSSM gauge-singlet scalar generates an effective 
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II parameter. For these reasons, having a nonvanishing /i' in the soft-breaking sector both 
solves the fi problem and serves as if there is a ;U parameter in the superpotential. 

The present work is organized as follows. In Appendix A we give the full list of renor- 
malization group equations (RGEs) for all rigid and soft parameters of the theory (as we 
hereafter call 'non-holomorphic MSSM' or NHSSM for short). In Appendix B we list down 
solutions of the RGEs of all model parameters as a function of their boundary values taken 
at the scale of gauge coupling unification Mqut ~ 10^^ GeV. An important parameter of the 
theory is the ratio of the Higgs vacuum expectation values: tan/? = {H^)/{H^). In solving 
the RGEs we will consider low (tan (3 = 5) and high (tan (3 = 50) values of tan /3 separately. 
In Sec. 2 we analyze the Z boson mass, in particular, its sensitivity to GUT-scale parame- 
ters. Here we will clarify the differences and similarities between the MSSM and NHSSM. 
In Sec. 3 we will discuss sfcrmion masses in the MSSM and NHSSM for the purpose of 
identifying their sensitivities to GUT-scale parameters, in particular, /iq and /Iq. Neutralinos 
and charginos are considered in the same section. In Sec. 4 we will discuss renormalization 
group invariants in the MSSM and NHSSM in a comparative manner so as to know what 
remains scale invariant in two distinct structures. In Sec. 5 we conclude the model. 

2 Fine-tuning oftheZ boson mass: MSSM vs. NHSSM 

It is well known that supcrsymmetry (SUSY) is not an exact symmerty of Nature, and there 
is no unique mechanism (gravity mediation, gauge mediation, anomaly mediation, etc.) for 
realizing its breakdown. From the viewpoint of Non-Standard Soft Breaking in the Minimal 
Supersymmetric Standard Model (NHSSM), on one hand, its predictions should reproduce 
the SM agreement with data, ensure unification of gauge couplings at the Grand Unified 
Theory (GUT) scale with minimal particle content, and on the other, it should preserve 
naturalness with soft terms [8]. 

It is expected that in the near future thanks to LHC and its successors, experiments 
related with superparticle masses and mixings will yield enough information to distinguish 
between various GUT-models and supcrsymmetry breaking mechanisms (see e.g. [9]). Taking 
gravity-mediation as the mechanism responsible for SUSY breaking, it is important to explore 
how the soft terms arc induced: holomorphic soft terms of the minimal model or those of 
the NHSSM with or without R parity violation [10]. In this work we will concentrate on 
NHSSM with exact R parity deferring the effects of R parity violation to a future work. 

Presently, apart from a number of observables in the fiavor-chaning neutral current sector, 
the Z boson mass is the main parameter that relates precision measurements to soft masses. 
In other words, the soft terms must self-organize so as to reproduce the measured value of 
the Z boson mass [8]. Hence, it is profitable to analyze Mz in the MSSM and NHSSM in a 
comparative fashion. 

2.1 Evolution of soft terms 

For the soft breaking parameters of the NHSSM [8] , we use one-loop Renormalization Group 
Equations (RGEs) [11] and thereby express their weak scale values in terms of GUT boundary 
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conditions (see Appendix A). Once weak scale mass values of SUSY particles are known, it 
will be possible to make educated guesses as to the GUT side. Meanwhile, the most general 
semi-analytic solution set of the RGEs for the NHSSM is too large for practical purposes 
to carry out phenomenological analyses which we present in Appendix B. Nevertheless, the 
number of free parameters can be considerably reduced if one assumes the universality of 
the soft terms at the GUT scale. In this case solutions are phenomenologically more viable 
and they can be found in Appendix C for all soft terms. Our choice for the GUT scale 
universality condition can be stated (dropping the contributions of all fermion generations 
but the third family) as some prototype structure inspired from minimal supergravity: 

At,,,M ^ A, , 4,,,(0) ^ A!, , Mi,2,3(0) ^ M . (4) 

Glearly, one may relax all or part of these conditions whereby obatining a larger parameter 
space augmenting the results presented in Appendix B. One should note that even if universal 
soft masses are assumed at the Planck scale, consideration of different boundary conditions 
for all soft terms including phases is more elegant, but then it gets difficult to achieve cetain 
clear-cut statements from the phenomenological side. To evade this cumbersome reality one 
needs certain inspirations which can be expected from string models. In order to use the most 
general one-loop solutions presented in this work, one can choose for instance, if the initial 
value of gauginos are not necessarily the same, then M30 7^ M20 7^ Mio can be implemented, 
and this approach can be generalized to all soft breaking terms. 

One of the most important distinctions is that, in the MSSM none of the soft masses 
depend on the initial value of ji, whereas in NHSSM both A! parameters and soft masses 
do depend on //q. Using the universality conditions of (4), let us present some of the soft 
masses in both of the models for low tanfi choice {tan/3— 5). In the MSSM masses of up 
and down Higgs at the weak scale can be expressed using boundary conditions of common 
gaugino mass, cubic and soft mass squared terms, 

m%Jtz) = -0.0S7 Al + 0.38AoM - 0.16ml - 2.8M^ 

m]j^{tz) = -O.OOSSAg + 0.011 AqM + 0.99m^ + 0.49M^ (5) 

whereas in the NHSSM also have primed trilinear couplings, 

m]j^{tz) = -0.087 Al + 0.1 A'^ - 0.16m^ - 2.8M^ + 0.067A[,^[, + O.Ufx'^ + O.SSAqM, (6) 
m]j^{tz) = -0.0033^^ _ o.37^(f + O.QQm^ + O.A9M'^ - 0.31^^ + + O.OII^qM. 

As it is seen in (5,6), at the electroweak scale, the results are the same except primed trilinear 
couplings and /iq, /Xq terms. As a matter of fact NHSSM predictions reduces to that of MSSM 
results under the following transformation: 

H', A[, 4, a; ^ , m]j^ ^ ^ m]j^ ^ + /x^ (7) 

which declares that NHSSM is a beautiful extension of the MSSM. In the NHSSM, notice 
that the contribution of Aq terms is not of the same order of Al terms for all soft masses. 
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hence trihnear and primed-trihnear couphngs are not symmetric (see Appendix C). What 
is more interesting is that, for both of the models, all soft masses depend heavily on the 
gaugino masses with the exception of leptons mf^ ^ . Among others mf^ is the most sensitive 
not only for gaugino masses but also for the initial value of /j,', for the latter mjj^ is the least 
sensitive in the NHSSM. 




[ - 
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t = (47r)-2 ln{Q/Qo) 

Figure 1: Scale dependency of gauginos in both of the models. Notice that here the boundary 
value of M is assumed to be 1 TcV. Scale dependency is expressed by dimensionless t such 
that to corresponds to 1.9 x 10^^ GeV . Here, Bino is at the bottom, followed by Wino and 
Gluino. Note that the same figure shows unification of gauge couplings. 

At this point it is appropriate to stress that there are also common model independent 
predictions like the evolution of gauiginos (i.e. see Fig.l), which stems from the insensitive- 
ness of gauge and Yukawa RGEs to both of the models at one-loop. On the other hand, 
trilinear couplings and other soft terms can be seen, in a way, to transformed into a new set 
in which /i terms are replaced with primed terms. 



2.2 Mz boundary 

For both of the models, as one of the most crucial constraints for the SM agreement with 
data, mass of the Z boson should be considered first, for a successful electroweak symmetry 
breaking. Notice that in the MSSM, in order to get the observed value of Mz, a delicate 
cancellation between the Higgs masses and /x is required, which is the famous n problem (see 
i.e. [12], [5]). Instead of ji parameter of the MSSM, NHSSM bears At',Ab',Ar' and ji' and its 
interesting effect can be seen by minimizing the scalar potential of the NHSSM which brings 
the constraint 

MUtz) ^ rnl^{tz)-tan^l3ml^{tz) 

2 tan2/5 - 1 ' ^ ' 

The Z boson mass depends on /ig rather strongly in the MSSM. As an example for tanj3=b, 
MSSM constraints can be expressed under the assumption of universality as 

^' = O.OM^ + 0.21mg + SM^ - 0.92//^ - 0.39^M. (9) 
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However, in NHSSM it does depend on ji' rather weakly e.g. a 10% change in //q generates 
only a 0.1% shift in M|/2. To make a comparison, in the NHSSM for the same value of 



MUtz) 



0.09A^ - 0.12^;; + 0.21m^ + 3M^ - 0.082AX - 0.12/x[f - 0.39/loM. (10) 



For the sake of visuahzation of the NHSSM and MSSM reactions we define dimensionless 
quantities ^i{tan(3) such that the Z constrain can be expressed as 



which can be used also for MSSM with obvious modifications. In the range tanl3 e [2,60], 
weights of 7's can be inferred from Figs. 2, 3 and 4. 




tan/3 tan/3 
Figure 2: Evolution of the coefficients of /x^ terms versus tanP e [2,60] in the MSSM (left), 
and of //'^ terms in the NHSSM (right) satisfying Mz constraint. 

In addition to relaxing sensitivity on the //q terms, we observe that tan(5 changes the 
sign of the ^' contribution in the NHSSM, and this situation has important consequences 
on the model building business. Note that in the MSSM contribution of terms is always 
destructive (assuming it is real), whereas by staring the oscillatory behaviour of //'^ with 
different choices of tan[3 (see Fig. 2) one can find a specific prediction for tan[3 such that /x'^ 
dependency of the M| completely vanishes in low and high regions, in addition to destructive 
or constructive contribution regions. Such special points can be called as turning points and 
this corresponds to ~ 49.25 for high tanfi in the NHSSM under the assumption of universal 
terms. Of course relaxing the universality assumption brings different turning points. 

In addition to capability of getting rid off n' terms for specific angles, NHSSM deserves 
new phenomenological approach which can be inferred from the Figs. 2, 3 and 4. In the figures 
tanf3 evolution of the coefficients of mass dimension 2 terms satisfying the Mz constraint is 
given. Behavior of other terms should also be taken into account, which can be performed 
using the Fig.4 to unfold the reaction of the model satisfying the mass boundary of the Z 
boson. To make a comparison with the MSSM we also presented the tan/3 evolution of similar 
coefficients in Fig. 3. Using figures presented here, one can find appropriate regions that 
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Figure 3: Evolution of the 71,3,4,7 terms versus tanP e [2,60] in the MSSM 

satisfy the Z constrain in both of the models for any tanf3 in the range [2,60]. Nevertheless, 
one should notice that constraints presented here are at the tree level and in the universal 
region, which might change when radiative corrections or anomaly boundary conditions are 
considered. 

Consequently, supersymmetry breaking with non-standard soft terms has an important 
virtue of reducing the sensitivity of M| to the initial value of the n parameter. However, in 
both cases, the MSSM and NHSSM, the Z boson mass exhibits a strong sensitivity of the 
gaugino masses. This follows mainly from the asymptotic freedom of color gauge group. 



3 Spectrum of sparticles in Minimal Supergravity: 
MSSM vs. NHSSM 

From the viewpoint of realistic model building approach any model should satisfy other 
collider bounds besides Mz, however we know from direct searches that no supersymmet- 
ric particle is observed yet, which can not set tight bounds on the spectrum of masses of 
SUSY particles [16]. Meanwhile mass of Higgs boson can be considered as on the verge of 
experimental verification if low scale supersymmetry really exists. We consider particle data 
group restrictions on the mass of sparticles and simply accept the lower bounds of LEP 2 
rusoft > 100 GeV, for the lightest chargino and neutralino half of Z boson width is accepted 
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Figure 4: Evolution of the 71^2,3,4,5,7 terms versus tan(3 e [2,60] in the NHSSM 



[17]. For simphcity and clarity, again, in this section we require all scalars to acquire a com- 
mon mass mo, all gauginos to be mass-degerate with M, all triscalar couplings to be and 
all non-holomorphic triscalars to be A'q all fixed at the GUT scale. In fact, suppression of the 
flavor- changing neutral currents as well as the absence of permanent electric dipolc moments 
already imply that the soft-breaking masses cannot be all independent and arbitrarily dis- 
tributed; they must be correlated by some organizing principle operating at the unification 
scale or above. With this assumption one can predict mass of lightest Higgs boson at tree 
level using the scalar Higgs potential of the NHSSM which brings the constraints 

mjj^ = mltan/S - (M|/2) cos2f3, (12) 
mjj^ = ml cot(3 + (M|/2) cos2/?. (13) 



During the numerical investigation, we look for real and positive soft terms in the range 
[0, 1000] GeV, which results in successful electroweak symmetry breaking patterns for low 
tanP option. In this case by noting the collider lower bounds on the mass spectrum, param- 
eter space can be restricted to a good extend, without additional assumptions (like no-scale 
[18], or some other string inspired models). With the same range proposed for GUT bound- 
aries there is no succesfuU candidate in high tan/S region, while the universality assumption 
of (4) in charge. When the electroweak symmetry is broken mass eigenstate of the lightest 
neutral scalar should satisfy > 114 GcV with radiative corrections. By expanding the 
scalar potential around the minimum trcc-lcvcl masses of the fields can be found as 

m^o = 2ml/sin2(3, (14) 

m]j± = m\o + M^, (15) 

ml\HO = ^{m\o +mlT + m|)2 - AMlm\oCos'^2(5), (16) 

when one-loop quantum corrections are considered SM like Higgs boson gets the largest 
contributions from t and b squarks. Notice that without quantum corrections mass of the 
lightest Higgs boson can not satisfy the experimental boundary, hence we study this issue in 
section 3.3 for NHSSM without GP violation; MSSM results including GP violation can be 
found in [19, 20] Analytic forms of mi^ and mf^ is given in the following subsection which 
will be needed in correction business. 



3.1 Sfermions 

For scalar fermions the relation between gauge eigenvalues and mass eigenvalues of the 
NHSSM particles can be read from the mass-squared matrices. Following that aim, we 
provide explicit expressions for the mass-squared matrices of squark and sleptons using 
reference [10]. The stop matrix is: 

/m|^+mf + i(4M2,-M|)cos2/3 mt{At - Af coi (3) \ 

V mt{At-At> cot 13) ml^ + mj-l{M^-M^z)cos2f3j ■ ^ 

for which we obtain the following eigenvalues 

= ^{6{2m^t+mtL + mln) + 3Ml;Cos2/3 (18) 

T sJaxcos2(5 (12(72 + (Ticos2(5) + 36 [4A|m| + cri + AAfmlcotp {-2At + A'tcotp)]}, 
where ai — 8M^ — 5Af| and (T2 — mfj^ — mfji. Similarly for the bottom squarks we have: 

f m^^ + ml -1{2M^ + M^) cos 2p mb{Ah - A^^ tan /3) \ , . 

V m,{A,-A,,tB,nP) ml^ + ml + 1{M^ - M^^) cos2p J ^ 

with eigenvalues 

= ^{6{2ml + m',^ + mlj,)-3Mlcos2p (20) 
T'^(73Cos2/3( 12(74 - (73COS2/3) + 36 [4:Alml + (j| + 4.Avmltanl3{-2A^ + A[tanl3)\}, 
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(21) 



where (73 = 4M^ — M| and (74 = m^^ — mf^. For the tau sleptons we have: 

il +ml-\ {2M^ - M|) cos 2/3 (A, - tan /3) 

rririAr - a; tan /3) mf^ + + (Af^ - M|) cos 2/3 

for which eigenvalues can be written as 

ml^^ = ^{2{2ml + ml + ml)-Mlcos2(3 (22) 

T\/cr5Cos2p{a5 - 4(76Cos2/3) + 4 [AA^ml + ai + AA'^mltan(3{-2Ar + A^ian/?)]}, 

where (75 = 4Af^ — 3M| and ue = — mfjij. Exphcit expressions related with each of the 
elements of these matrices can be extracted from the Appendix C of this work for low and 
high tanf3 choices. In the MSSM sfermion masses depend on yUo only via their (1,2) and (2,1) 
entires whereas in the NHSSM /Iq appears in all entires including (1,1) and (2,2). When all 
the Yukawa couplings are set to zero, except ht and /i^-, it is interesting to observe SUSY 
loop effects on the mass squared terms (see [26] and [27]). 



3.2 Charginos and Neutralinos 

The last step is to compare the mass eigenvalues of neutralinos and charginos. Neutralino 
values can be read from the following matrix, which resembles the mixing of Higgsinos and 
neutral gauginos 



/ Ml —Mz cosf3 sinOw Mz sin/3 sinOw \ 

M2 Mz cos(5 cosOw —Mz sin/3 cosOw 

—Mz cos/3 sinOw Mz cos(3 cosOy/ —/i' 

\ Mz sinj3 sinOw —Mz sin(3 cosOw -/J'' / 



(23) 



Similarly charginos are mixtures of charged Higgsinos and charged gauginos with the mass 

M2 V2Mwsinp\ , . 

V2Mwsin/3 /i' J ' ^ ^ 

Since we assume R-parity conservation LSP is the lightest neutralino. Explicit form of matrix 
elements can be found in Appendices for low and high values of tanp. 



3.3 Higgs boson mass and LEP bounds 

In this section we will compute the Higgs boson mass in NHMSSM. The main impact of 
the non-holomorphic soft terms on the Higgs boson masses stems from the modifications 
in the sfermion mass matrices. Indeed, as one infers from the forms of the sfermion mass- 
squared matrices in Sec. 3.2, the mixing between the left and right-handed sfermions are 
described by the holomorphic triscalar coupling A^ and the non-holomorphic contribution 
A'^. The left-right mixing thus changes from flavor to flavor in contrast to MSSM where A'j: 
is replaced by flavor-insensitive quantity /i parameter. 
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For a proper understanding of the Higgs sector it is necessary to implement the loop 

corrections as otherwise the tree level masses turn out to be too low to saturate the experi- 
mental bounds. The radiative corrections to Higgs boson masses and couplings have already 
been computed in [19, 20] including the CP-violating effects. Concerning the neutral Higgs 
sector, it is useful to use the parametrization 



(01 + i(pi) ; = (02 + i(p2) , 



(25) 



where 0i_2 and ipi^2 are real fields. The Higgs potential, including the Coleman- Weinberg 
contribution, reads as 



Higgs — 



-m 



2 /2 

t 9 + 9 (\ uQ\2 |o-0|2\ 



1 



+ 



647r2 



Str 



(26) 



where g and g' stand for the SU(2) and U{1)y gauge couplings, respectively [g'"^ = fs'i). Qo 
in (26) is the renormalization scale, and A4 is the field-dependent mass matrix of all modes 
that couple to the Higgs bosons. The masses of the quarks are to be taken into consideration 
of which the most important contributions come from: 



nil. 



l>4 



2 I ,o2 
1 + 9^1 



my. 



2 * 



(27) 



Now, using the eigenvalues of the field-dependent squark mass matrices (18,20) in (26) one 
can systematically compute the Higgs boson masses at the minimum of the potential obtained 
via the conditions 



dV, 



Higgs 



0(1)1 



= 0, 



dV, 



Higgs 



= 



with (ifi) — {if 2) — and 



^1)' + 



Ml 
9' 



~ (246 GeV)^ 



(M 



tan/5. 



(28) 



(29) 



where g'^ — {g'^ + g'"^) / A. The mass matrix of the neutral Higgs bosons are computed from the 
matrix of second derivatives of the potential (26). Notice that after including the one-loop 
corrections to the Higgs potential, the Z mass becomes dependent on the top- and stop quark 
masses too [29]. In this case there will be a correction term 



MUtz) _ mjj^(tz)-tan''Pmjj^(tz) 



A%{t,b) 



where 



tan^P - 1 



f{ml)-f(ml) 

+ 2m^ + /(mjj + f[m^^, 



m 



t2 



(30) 



(31) 
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and 



(32) 



Similarly A|(6) can be found with the t ^ h substitution. This corrections require a large 
amount of fine tuning if the mass splitting between the particles and sparticles is large [8]. 

The Goldstone boson — (pi cos P — (p2 sin /3 is swallowed by the Z boson. We are then 
left with a squared mass matrix Ai^ for the three states (p — (pisinP + (p2 cos /3, 0i and 02- 
If the theory has CP-violating phases (via the phases of the triscalar couplings and n') the 
(/? mixes with 0i and 02. In the CP-conserving limit, however, ip decouples from the rest, 
and assumes the mass-squared: 



Ml 



_ 2 _ 2m| 2 
aa-"^^" sin (2/3) ^ sin(2/3) 



where 



F{mi,m2) = 



hfAtAt,F{ml,ml) + hlAbAvF{m\,ml^) 
3 /(m?)-/(mi) 



327r2 



2 2 
777.2 ~ ^1 



(33) 



(34) 



The remaining real scalars 0i and 02 mix with each other via the mass-sqaured matrix: 



H 



+ 



+ 



M| cos^ l3 + m\ sin^ 

g{m\,m\)Rt [hlRt - coi(5X^ + cot/3i?t log 



3m2 
8^ 

3mg 
8^ 



2 2 

mt 



1 2 
01 62 



2 

I2_ 

2 

ti 



(35) 



+5(m?^, m?Ji?^ (/i^i?^ + X,) + log 



= M| sin^ /3 cos^ /3 



+ 



6i 



3m2 
"8^ 



2 2 

mf mf 



Qt 



2 

+9{ml,ml)R', {h',R', + X,) + log ^ [x, + hh^ - g') R', 

tl 



Xb + {2}i, 




3ml 
8^ 



g{ml^ , ml^)Rb (hlRb - ternpXb) + tan/3i?6 log 



2 

mf 



m 



bi 



(36) 



where 



and 



X, 



/ 2 2\ r. + 1 "^1 

^(mi, 7^2) = 2 2 2 log 



2 2 

TTlf — 777,2 



- 3g^ 
12 



2 2 
?7i~ — mf 

12 tl 



,/2 



2 ' 
777^ 



3^^ 



12 



777- 
02 



777- 
01 



(37) 



(38) 
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Rh 



Al cot/3 + AtAf 



R' 



Af + AtAf cotp 



t2 



m 



ti 



m 



t2 



m 



ti 



_ Al tan/j + Ai,Ay , _ Af + AfoAy tan/3 



ml 

02 



ml 

01 



2 

mi 

02 



2 

mt 

01 



(39) 
(40) 



It is known that, the two- loop corrections to Higgs boson mass are reduced at the renor- 
mahzation scale Qq — mt hence our choice hereon. To give a concrete example of NHSSM 




Figure 5: Scale dependence of the couphngs ht (left) and of hi, (right) for different choices 
of tan(3 e [5,55]. In both of the figures topmost curves correspond to tan(5—hh. 

benchmark we now list mass predictions of the model for low tanfi with the input param- 
eters; fht{tz) = 170, fhi,{tz) = 2.92 and m^(tz) = 1.777 GeV and take the GUT boundary 
values of soft terms as the following set 

M = 160, mo = 683, = 400, Ao = 800, A'^ = 1000, m^o = 430 (41) 

which brings the following predictions 

mi^{tz)= 291, m,-^(tz) = 626, 
m^^(tz) = 600, m^^(tz) = 791, 
^Ti{tz) = 683, mr2{tz) = 695, 

"^xU3,4M=63, 120, 392, 407, (42) 

m^±^{tz) = 119,407, 

mAo{tz) = 289, mH±{tz) = 300, 

rriHoitz) = 291, muoitz) corrected = 123, 

where all masses are given in GeV. 

Since NHSSM covers MSSM any prediction of the classical MSSM results can be repro- 
duced in non-holomorphic case with the appropriate boundaries. But the extension enriches 
us with more opportunities. What it is important here is the degree of freedom offered by 
NHSSM. As it was stressed in [10] for mo > M it turns out that < 0.4 M in the MSSM 
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whereas in the NHSSM this constrained is significantly relaxed. Note that in our example we 
assumed all soft terms as if they are real and positive without considering any specific model, 
whereas one can study i.e Aq = —M which arises in certain string inspired models. Under 
the light of these observations, it should be stated that, NH extension of the MSSM not only 
covers the classical MSSM but also offers novel features that can ease the shortcomings of 
the MSSM, which should be studied in more detail. Actually, in addition to LEP limits on 
the SUSY mass spectrum, one should also deal with the constraints from b 57 decay (as 
we do in next subsection) and the lower limit on the lifetime of the universe, which requires 
the dark matter density from the LSP not to close the universe on itself [30] . 



3.4 6 — ^ S7 Decay 

Presently, one of the most accurate observables which can severely constrain the soft masses 
is the branching ratio for the rare radiative inclusive B meson decay, B — > Xs^. The main 
interest in this decay drives from the genuine perturbative nature of the problem and also 
from the striking agreement between the experiment and the SM prediction. Indeed, the 
measurements of the branching ratio at CLEO, ALEPH and BELLE gave the combined 



13 



result [21] 



BR {B ^ X,7) = (3.11 ± 0.42 ± 0.21) x lO"'' (43) 

whose agreement with the next-to-leading order (NLO) standard model (SM) prediction 
[22] 

BR (B ^s7)sM = (3-29 ± 0.33) x 10"^ (44) 

is manifest though the inclusion of the nonperturbative effects can modify the result slightly 
[23]. That the experimental result (43) and the SM prediction (44) are in good agreement 
shows that the "new physics" should lie well above the electroweak scale unless certain 
cancellations occur. 

The branching ratio for B Xg'y has been computed up to NLO precision in the 
MSSM [24]. The W boson and charged Higgs contributions are of the same sign and thus 
the chargino-stop loop is expected to moderate the branching ratio so as to respect the 
experimental bounds. The recent measurements of BR (5 ^ Xsi~^i~) [25] imply that the 
sign of the total b s'-f amplitude must be same as in the SM. This eliminates part of the 
supersymmetric parameter space in which the total amplitude approximately equals negative 
of the SM prediction. In spite of these, however, the present experimental results do not 
exclude stop masses around a few Mz as long as At and A!f. are of opposite sign [24] . 




-0.2 -0.15 -0.1 -0.05 -0.2 -0.15 -0.1 -0.05 



t t 

Figure 7: A sample plot of the scale dependence of the trilinear couplings for tanf3—5 (left), 
tan(3=50 (right) with Aq = A'q=0, M=150 GeV and //'=1000 GeV, which show a candidate 
region where At and A'^ are of opposite sign. 



To accommodate differing signs of trilinear couplings in the NHSSM we present another 
example using the following input parameter 

M = 200, mo = 787, h'q = 400, Aq = 900, A'q = -1500, mgo = 414 (45) 

which yields the following predictions 

mt^{tz)=362, mi^{tz) = 728, 
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rrii^itz)^ 711, mg^(iz) = 930, 
mr,{tz)^ 787, mr^{tz)^801, 

m^o^^^ (tz) = 79, 150, 392, 409 (46) 

m^±^{tz) = 149,408, 

mAo{tz) = 299, mH±{tz) = 310, 

mHo{tz) = 301, mho{tz) corrected = 120, 

here again all masses are given in GeV. 

4 Renormalization Group Invariants in the MSSM and 
NHMSSM: A Comparative Analysis 

Renormalization Group Invariants (RGIs), which can be used to relate measurements at the 
electroweak scale to physics at ultra high energies provide important information about high 
scale physics due to the scale invariance of the quantities under concern [31, 32]. Since the 
coupled nature of the RGEs disturbs analytical solutions it would be beneficial to know if one 
can construct certain invariants that give relations among the spectrum of supersymmetric 
particles. Indeed, RG invariants may provide a direct, accurate way of testing the internal 
consistency of the model and determine the mechanism which breaks the supersymmetry. 
Such quantities prove highly useful not only for projecting the experimental data to high 
energies but also for deriving certain sum rules which enable fast consistency checks of the 
model. Assume there is a measurement which tells a specific relation between some of the 
soft masses, then, it can be easily probed whether this relation survives at different scales or 
not, with the help of scale independent relations, which in turn shows the way how SUSY is 
broken. 

In this part we will discuss RG invariant observables in supersymmetry with non-holomorphic 
soft terms and compare with existing MSSM results with the assumption that there is no 
flavor mixing and soft terms obey the universality condition mentioned previosly. Neverthe- 
less, it should be kept in mind that we study one-loop RGIs which differs when R parity or 
higher loop effects are taken into account. 

To begin with, note that lagrangian of the NHSSM (2) has parameters deflned at a specific 
mass scale Q which can physically range from the electroweak scale Q = Mz (the IR end) 
up to some high energy scale Q = Qo (the UV end). For determining the scale dependencies 
of the parameters the RGEs are to be solved with proposed boundary conditions either 
at IR or UV. In what follows we will write them in terms of the dimensionless variable 
t = (47r)^^ ln(Q/(5o), and solve for the parameters in terms of their UV scale values by 
taking into account the fact that the gauge and Yukawa (at a given tan/?) couplings are 
already known at IR end. 

We should deal with the rigid parameters in both of the models as a first step. The RGEs 
for gauge and Yukawa couplings form a coupled set of first order differential equations and 
can be found elsewhere (i.e. see [11]). Now one can solve them at any scale at one loop order 
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without resorting to other model parameters. However, expanding this set of equations by 
including the RGE of the jj! parameter one finds that 



is a one-loop RG-invariant. For the classical MSSM invariant fi' ^ fi substitution suffices 
(MSSM was also mentioned in [31]). Here the powers of the Yukawa and gauge couplings 
follow from group-thcorctic factors appearing in their RGEs. This invariant provides an 
explicit solution for the n' parameter 

27 21 10 256 _9_ 73 

vrt 'im ( ''<W\" ('>i(t)\"' (i'r{t)Y- (giio)"]'" (gimV" /siWV™ usi 

once the scale dependencies of gauge and Yukawa couplings are known either via direct 
integration or via approximate solutions the RGE of the /x' parameter involves only the 
Yukawa couplings, g2 and gi though this explicit solution bears an explicit dependence on 
g^. This follows from the RGEs of the Yukawa couplings. One of the most interesting sides 
of this invariant is that weights of all gauge and Yukawa couplings is made obvious. With 
this equation one can determine the amount of fine tuning to satisfy Z mass boundary (see 
reference citedurmus for a detailed discussion on this issue). Another by-product of the 
invariant Ii is that the phase of the fi parameter is an RG invariant. Since the contribution 
of higher order loop effects affect invariance relation of (47) ~ 2 — 3%; an effect likely to 
get embodied in the experimental errors encourages us to work at one-loop order. On the 
other hand, once the fiavor mixings in Yukawa matrices are switched on there is no obvious 
invariant like (47) even at one loop order. 

We continue our analysis with the construction of the RG invariants of the soft parameters 
of the theory. Of this sector, a well-known RG invariant is the ratio of the gaugino masses 
to fine structure constants 



with one-loop accuracy. This very invariant guarantees that 



so that knowing two of the gaugino masses at Q = Mz suffices to know the third - an 
important aspect to check directly the minimality of the gauge structure using the exper- 
imental data. Related with this invariant it is useful to state the well known mass ratios 
M^{tz)/M2{tz) = 3.46 and M2{t z) / Mi{t z) = 1.99 at one-loop order. The invariant (49) 
pertains solely to the gauge sector of the theory; it is completely immune to non-gauge 
parameters. At two loops I2 is no longer an invariant; it is determined by a linear combi- 
nation of gaugino masses and trilinear couplings. Combining (48) and (50) one concludes 




(47) 




(50) 
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that the chargino and neutrahno sectors of the theory are connected to the UV scale via 
the gauge and Yukawa couphngs alone. The equation (50) suggests that M'^{tz)/M'^{0) is 
much larger Mi^2(^z)/^i,2(0) due to asymptotic freedom, and these coefficients stand still 
whatever happens in the sfermion and Higgs sectors of the theory 

A by-product of the invariant (49) is that the phases of the gaugino masses are RG 
invariants (like that of the fj, parameter). However, this is correct only at one-loop level; 
at two loops the phases of the trilinear couphngs disturb the relation between IR and UV 
phases of the gaugino masses. 

Another invariant of mass dim-1 is related with the B parameter for which we obtain: 

„ 27 , 21 10 , 256 _ 9 _ 73 _ 

U ^ B At Afe M3 M2 H Ml 

^ 61 61 61 183 61 2013 

+ CiA; + + CsA; - (Ci + C2 + C3)/x', (51) 

with arbitrary coefficients Cj such that in the limit ^ ^, /x' — > /i it reproduces the well known 
MSSM invariant which can be expressed in terms of other parameters 

27 21 10 

Bit) = B{0) + -iAit)^AM) + ^iMt)-MO)) + ^iA^{t)-AM) (52) 

+ -M3(0) f 4^ - lU ^Af.(O) f 4^ - ll - ^M,(0) f - iV 

183 ^ \g,{Oy J 61 ^ [g^iOy J 2013 ' \g,{Qf ) 

Concerning mass dimension-2 terms we obtain a general invariant relation in the NHSSM 
by brute force as follows 

Cl , 9C2 , C3 ^ Ci\ 2 . f—Ci ^ 3C2 C3 



^ _ ^ _ -L 2 , 9C2 C3 3C4 



+ - "Jill - + ':2:i\rnUt) + - — - - - — ]mi it) 

\ 2 16 2 2 ; ^ V 2 16 2 2 >' ^ 



33 44 



+ c,A'^{t) + c,A'^{t) + (i) _ + C5 + C6 + C7) /x'^(t). (53) 

where q are arbitrary constants. To visualize our results lets set all coefficient to zero but 
C5,6,7 we then obtain 

cs^f (t) + c^A'^it) + c^A'^it) - (C5 + C6 + c^)ii'\t), (54) 

which is obviously invariant in the limit A'^^^, fi' ^ jjL. Note that using this limiting case one 
can obtain another invariant, when supplemented with m^^^lt) — > m\^^{t) + ^x'^{t) brings 
the most general form of MSSM invariant mass of dim-2. In the cases when we relax these 
substitutions we obtain more general structures. Now we vary the coefficients of various soft 
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masses for constructing invariants in terms of Mj and // parameters. Using this freedom, 
when we set Ci = —3, C4 — I and all other coefficients to zero, we get 

h = -2ml{t) + mUt) + 3 \M2{t)f - ^ \M,{t)f (55) 
and similarly various patterns of the coefficients give rise to 

h = ml,^(t)-^m^,(t)-<(t) + ^|M3(t)r-^|Mi(t)|^-|/.'(t)f, 

^8 = ml^{t)+ml^{t)-2ml{t)-m2{t)? + ^\M,{t)\\ (56) 
^9 = m^jt)+m^^(t)-3m2jt)-mfjt) 

+ ^ |M3(t)|' - 3 \M,{t)\^ + 1 \M,{t)\^ - 2 \^x\t)\' , 

+ \ \M,{t)f - \ \M,{t)\' + ^ |Mi(i)|^ - \^,'{t)f , 

which should be compared with the results of (see [32]). Clearly, one can construct new 
invariants by combining the ones presented here or by varying the coefficients expressed as 
Cj. Although the results presented here and the results of [32] coincide a term is observed 
to be missing in some of the invariant equations. This stems from the definitions and 
frameworks i.e. we work within minimal supergravity (with non-holomorphic soft terms). 
Here we confirm the results of [31, 32] in certain limits and we also generate new invariants. 

The general form (53) and the invariants that follow could be very useful for sparticle 
spectroscopy [16] in that they provide scale-invariant correlations among various sparticle 
masses. All the invariants presented here show non-anomalous behaviors unless they bear 
n' terms. As an example lets take Jg Fig. (8), which demonstrates the fixed behavior. No- 
tice that while it is scale-dependent, it is still very useful since its dependency is very soft. 
However, notice that they are obtained without noting flavor mixing and in the mSUGRA 
framework. Nevertheless, using them one can (i) test the internal consistency of the model 
while fitting to the experimental data; {ii) rehabilitate poorly known parameters supplement- 
ing the well-measured ones; {in) determine what kind of supersymmetry breaking mechanism 
is realized in Nature; and finally {iv) separately examine the UV scale configurations of the 
trilinear couplings as they do not explicitly contribute to the invariants. 

Consequently, if one single invariant is measured then all are done, and in case the experi- 
mental data prefer a certain correlation pattern among the invariants then the corresponding 
UV scale model is preferred. In this sense, rendering unnecessary the RG running of indi- 
vidual sparticle masses up to the messenger scale, the invariants speed up the determination 
of what kind of supersymmetry breaking mechanism is realized in Nature. 
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t = (47r)-2 ln{Q/Qo) 

Figure 8: Fixed point behavior of the anomal Ig against scale. Here we assume same weight 
for all soft terms (~ 40 GeV) and re-scale the figure (initial value of this invariant is ~ 
32TeV'^). 

5 Conclusion 

It is important to explore the features of MSSM and its extensions as general as possible. This 
will be clear as experimental data accumulates about the masses of all predicted particles, 
and for the time being it should be calculated at low energies using the RGEs. For that aim 
NHSSM offers novel opportunities which should be studied in more detail. Compared with 
its enrichments, there are not enough papers in the literature about the phenomenological 
consequences of the NHSSM. So we try to cover this issue from many sides. Because we do 
not know the mechanism of supcrsymmetry breaking, extensions of the MSSM should be 
taken seriously to ease the shortcomings of the MSSM. In this paper we explored the main 
features of NHSSM with minimal particle content and observe that, in addition to mimic 
the reactions of the MSSM (like gauginos or Yukawa couplings), NHSSM offers interesting 
opportunities. Even, under certain assumptions, it is possible to completely get rid of famous 
fi problem in the NHSSM, and this corresponds to two special turning points in low and high 
tan(3 regimes, which is not possible in classical MSSM. The price that must be paid is, facing 
additional primed trilincar coupling and fine tuning of parameters for GUT boundaries. 

One of the main results of this work is to present semi-analytic solutions of RGEs of 
NHSSM which enables one to study the phenomenology in detail. Using the solutions pre- 
sented here one can investigate the reaction of the NHSSM deeper. Notice that the solutions 
presented in the Appendices have nonzero phases which should be used to go deeper in the 
phenomenology. 

Another result is to present a general form of RGls which can be used to derive new 
relations in addition to those existing in the literature. We observed that by using existing 
RGEs one can construct RGIs with a simple computer code which indeed offers a very 
practical way of handling the equations. These invariants turn out to be highly useful in 
making otherwise indirect relations among the parameters manifest. Moreover, they serve as 
efficient tools for performing fast consistency checks for deriving poorly known parameters 
from known ones in course of fitting the model to experimental data, and for probing the 
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mechanism that breaks the supersymmetry. 
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A Explicit form of RGEs of the NHSSM 

For the NHSSM one-loop renormalization group equations can be found in [10] we also 
present here for the sake of completeness. 



+ 6hl Al - SChIJ^'^ - QglMl - 2g'^M^, (57) 



(3^2^^ = 6hl{mjj^+A^,+ml+ml) + 2hlAl, + 6hlAl, 

- 8CHfJ^'^-GgiMl-2g'^Ml (58) 

Pmi = (hi + ^hl + 3h^t)ml + 2hl A^,A^ + QhlAyA^ + Qh] At, At 

- 4:CHml + 6glnfM2 + 2 g'^Mifj,', (59) 

P^, = {hl + 3hl + 3h^t-'^CH)fJ>', (60) 

/3a,, = {hl-3hl + 3h^t)Ar' + QhlA,, + {AA,,-8ix')CH, (61) 



/3a. 


= 8hlA^ + QhlA^ + QglM2 + Q g'^Mi, 




(62) 


Pa, 


= {-hi + 3hl + hl)Ab, + 2Ar,hl - 2hl{At> - 2/x') + {AA^, 


- 8/x')Cj^, 


(63) 


Pa, 


32 14 
= 2hlAr + 12hlA^ + 2hlAt + —glM^ + QglM^ + y^^Mi, 


(64) 


Pa, 


= {hi + hl + 3h^t)At' - '2Ayhl + Aii'hl + (4A' - ^I^')Ch, 




(65) 


Pa, 


32 26 
= 2/^2^6 + l2hlAt + y^^Ma + 6 glM^ + -g-y"Mi, 




(66) 


^< 


= 2hl{ml + m^^ + m^^ + Al + Al - 2fi'') 
+ 2/1? {ml + ml + m\j^ + A?, + A? - 2/x'2) 
- ^^glMl-^glMl-lg^Ml 




(67) 


/3m2 

*R 


= 4/i?(m?^ + ml + m^^ + A?, + A] - 2//'^) - ^glMj - 




(68) 


Pml 


= 4/i,2(m,\ + ml + m^^ + A^, + Al - 2/i'2) - ^^fM^ - 




(69) 
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(3^.^ = 2hl{ml-rml-rml^-rAl,-rAl-2^,'')-&glMl-2g''Ml (70) 
= Ahl{ml+ml + m\,^ + Al,+Al-2^,'')-%g''Ml (71) 

/3m, = '2hM,gl (72) 

here 61,2,3 = (f , 1, -3), 5'' = |^?, = + l^gl Mgut = 1-4 x lO^^ GeV and Mz < 
Q < Mgut- By assuming that the SUSY is broken with non-standard soft terms; we 
obtained semi-analytic solutions for all soft terms through the one-loop RGEs given above 
and express our results at the electro-weak scale in terms of GUT scale parameters. Our 
results are presented for moderate {tan(3=5) and large {tan(3=50) choices. 



B Solutions of mass squared & trilinear terms in the 
NHSSM 

Using low {tanP = 5) and high (tan[3 = 50) values of tanP, the most general form of 
the mass-squared and trilinear terms can be written in terms of boundary conditions of 
gauge coupling unification scale which is roughly Mqut ^ 10^^ GeV. Notice that our phase 
convention is to assign 1, 2, 3 and 4 for Mi, M2, M3 and for other quantities it is obvious 
and can be inferred from the multipliers. 



B.l Low tanP regime 

m]j^{tz) = 0.000216A^^ - 1.59 x lO-^^^ A^^cos^b^ - O.OOOO2O3AfeoMioCos0bi 



- O.OOO19lA6oM2ocos062 - 0.000857AoM3ocos(/)63 - 0.00124A^, 

+ 1.73 X IQ-^ Ay^Ar^COS(t)b'r' - 0.000563^;,/ //[,COS0fe/4 - 0.0869^^ 




2 

to 

8 



+ O.OOOO648y4t„y46„cos0tb - 2.05 x IQ'^ At^Ar^coscjitr + O.OlO9AjoMiocos0a 

+ O.O672^„M2ocos0i2 + QM2At,M^^cos(t)u - 7-96 x 10"^^^ 

+ 2.25 X IQ-^ Ar^MiQC0S(pri + 1.19 x lO~^A^oM2oCos0^2 + 4.14 x IQ''' A^^M^ocos(prz 

- 0.000287^^, - O.OOO248yl^^/ioCos0^/,4 + 0.105^4^, 

- 0.000284At' Ay cos<i)t'w + 2.25 x lO^Mf A^>cos(l)t>r' + 0.0674^^/ HQC0S(l)t'4 
+ O.OOlOGMi^o - O.OO58MioM2ocos0i2 - O.O291MioM3ocos0i3 

+ 0.187M2o - O.2O6M2oM3oCos023 - 2.79M|o 

-f- 0.000217m^^o + 0.000217m^^o + 0.612m^^o 

- 7.98 X 10"^mf^o - §• x 10"^rnf^o - 0.388rr4o 

- 0.388m^^o + 0.136/x[f , (73) 
mH^{tz) = -0.0032^^^ + 5 X 10~%o A^cos^b^ + O.OOO18A6oMiocos0;,i 
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+ Om22Ab,M2ocos(j)b2 + OMAb^Msocos(j)b3 + 2.9 x 10"^^^^ 

- 2.4 X 10~^ Ay^Ari^cos(f)b'r' - O.OOO17^6^//()Cos06/4 + 0.00008^^^ 

+ O.OOO58Ao^6oCos0t6 - 4.7 x IO^^AqA^oCos^j^ - O.OOOO28Ao-^iocos0ti 

- .00029 At,M2ocos(f)t2 - 0. 0013 At,Msocos(j)t3 - 0.00078^1^ 

+ 0. 00018 AroMioCOS(j)rl + 0.0005A^M20COS(t)r2 - 7.9 X 10~^AroM3oCOS(l)r3 

+ 5.2 X 10-^^^, + 3.5 X 10"%^//oCos(/)^/4 - 0.37A| 

- 0.00026At>^Ah'^cos(l)t'b' + 7.5 x 10~^Aj, ^r^cos^t/^/ - O.31A[,AioCos0t/4 
+ 0.037Mi^o - 0.00013MioM2oCOS(/)i2 - O.OOO3MioM3oCos0i3 

+ 0.48M|o - O.OO4M2oM3ocos023 - 0.026M|o 

- 0.0032mg^o + ^i,o + 0.00029ml,„o 

- 0.00079mf^o - 0.00079mf^o - 0.0029mt\o 

+ 0.00029mj2^o + 0-6Afo > (74) 

= -0.00099^1^^ + 9.8 X 10"'^ AboAroCos(f)bT + 0. 000053 AboMiocos(f)bi 

+ 0. 00068 AboM2ocos(f)b2 + O.OO3A6oM3ocos063 - 0.00041A^j^ 

+ 3.1 X 10~^ Ab'^A^^cos(f)b'r' - 0.0002AAb'yQCos(f)b'4 - 0.029A^^ 

+ 0. 00022 AtoAboCos(j)tb - 1.2 x 10''^ AtoAroCos(j)tr + O.OO36Ao-^iocos0ti 

+ 0.022At,M2ocos(l)t2 + 0.1 At, M3ocos(f)t3 + 4.9 x lO'^A^ 

- 1.3 X 10~'^AroMioCOS(l)rl - 6.4 X lO~^A^oM2oCOS0^2 - 1-9 X 10"*^ AroM3oCOS(l)r3 

+ 0.000039^1^, + O.OOOO24A^///()Cos0^/4 - 0.089^^^^ 

- 0. 00018 At'^Ab'^cos(l)t'b' + 7.7 x lO^^^A^^cos^f^' - 0. 08 At'^i^'ocos(j)t' 4 

- 0.0081Mi2o - O.OO2MioM2ocos0i2 - O.OO98MioM3oCos0i3 
+ 0.38M|o - 0.07M2oM3ocos(/)23 + 5.4M|o 

- 0.00099m^^o - 0.00099m|^o - O-lSm^^o 

+ 4.9 X 10"^mf^o + 4.9 x 10"^mf^o + 0.87m^^o 

- 0.13m?^o + 0.3//(f, (75) 

= 0.00014^^^ - 1.1 X 10-'^ Abo AroCos(f)bT - O.OOOOUAboMwcos(f)bi 

- 0.00013A6oM2ocos(/)62 - O.OOO57^6oM3ocos063 + 0.00037A^j^ 

- 3.2 X 10-^Afe/ A^/cos(/)6v + 0.000042 A6/^//()Cos06/4 - 0.058A^^ 

+ O.OOOO43Ao^6oCos0t6 - 1-4 x 10"^^^ A^^cos^t^ + 0. 0072 AtoMiocos(l)n 

-8 a2 



+ O.O45AtoM2ocos0t2 + 0.2At,M3ocos(l)t3 - 5.3 x lO^M 

+ 1.5 X 10~^A^Miocos(l)ri + 8.X lO~^A^oM2oCos0^2 + 2.8 x lO~^A^(,M3oCos0^3 

+ 0.000077^^, + O.OOOO48A^^/XoCos0^/4 - 0.18Al^ 

- 0.000073A^,A^ cos0t/6/ + 7.7 x 10~^At>^Ar^cos(j)t'r' - O.16Aj,A«oCos0t'4 
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+ OM3M^o - O.OO39MioM2ocos0i2 - 0.019MioM3ocos(/)i3 

- 0.2M|o - O.14M2oM3oCos023 + 4.4M|o 
+ 0.00014m^^o + 0.00014m^^o - 0.26m^^o 

- 5.3 X lO'^mf^o - 5.3 x lO'^mf^Q - 0.26mt\o 

+ 0.74m? . + 0.6//(f, (76) 



ml^itz) = -0.0021A?^ + 2.1 x IQ-^Abo^^cos^fc^ + O.OOO12^oMiocos06i 

+ 0.0015^oM2ocos(/)62 + .0066 AboM3ocos(f)b3 - 0.0012^?^ 

+ 9.5 X 10-'^ Ab'^Ari^cos(f)b'r' - O.OOO53^6j^//oCos06/4 + 0.000053^^^^ 

+ 0. 00039 AtoAb^cos(l)tb - 2.2 x lO~%oA^oCos0t^ - 0.000019AtoMiocos(/)ti 

- O.OOO19Ao^2ocos0t2 - O.OOO86Ao^3ocos0t3 + 1 x lO'^A^^ 

- 2.7 X 10"'^ AroMioCOS(l)rl - 1.4 X 10~'^AroM2oCOS(l)r2 " 4.1 X lO~^A^oM3oCOS0r3 

- 4.5 X lO~^Al^ + 2.3 X 10~^ Ar;^n'ocos(l)r'4 + 0.00068A?^ 

- 0. 00029 At'^Ah'^cos(l)t'b' + 1.5 x lO~^A^,^r^cos0t/^/ + O.OOO38AiA*oCos0t'4 
+ 0.01771^1^0 - 0.000042MioM2oCOS(/)i2 - O.OOO2MioM3oCos0i3 

- 0.0017M|o - O.OO27M2oM3ocos023 + 6.3M|o 
+ lm?^o- 0.0021ml,^o + 0.0002m^^o 

+ 1 X 10-^mf^o + 1 X 10~^w?«o - 0.0019mt\o 

+ 0.0002m?^o + 0.0029//(f , (77) 



m]^{tz) = 9.6 X 10~^A^^ + 1.9 x lO^^AbgA^oCos^b^ - 3.1 x lO^^gMiocos^bi 

- 1.3 X lO~^^oM2oCOS062 - 1.8 X lO~%oM3oCOS063 - 1-3 X 10'^ Al^ 

+ 7.9 X 10'"^ Ab'^Aricos(f)b'r' + 5.4 x lO~^^bj^//oCos06/4 - 2.6 x 10"^^^^ 

- 1.3 X 10-"^ At, Ab^coscptb - 1.3 X 10''^ At, ArQCos(f)tr + 2.5 x 10~^ At,Miocos(l)n 
+ 1.1 X lO-%oM2ocos0t2 + 2.1 X 10-^Ao^3oWt3 - 0.00079A^ 

+ O.OOO18^^oMioCos0^i + O.OOO5>l^M2oCos0^2 - 1-8 x lO~^A^M3ocos0^3 

- 0.0004^^, - O.OOO32A^^//[)Cos0^/4 + 0.00018^^^^ 

+ 6.9 X 10~^ At'^Ab'^cos(l)t'b' + 6.8 x 10"^ Aj,^T^cos0tv + O.OOO1A^,/Wocos0t'4 

+ 0.038Mi2o - O.OOOO66MioM2ocos0i2 + 3.2 x 10'^ MioMsocoscpn 

+ 0.48M|o + 1.5 X lO-^M2oM3oCos023 + 4. x 10-^M|o 

+ 9.7 X lO^^m^^o - 0.00079m|^o - 6-6 x 10~^m^^o 

+ Imf^o - 0.00079mf^o + 9- x 10-^,0 

- 6.6 X 10-^m?^o + 0.0012//(f , (78) 
mf^{tz) = 1.9 X lO^^A^o + 3-9 x 10~^ Ab,AroCos(l)br - 6.3 x lO~%oMiocos06i 



24 



- 2.6 X lO-^AoM2ocos062 - 3.7 x 10-^AboM3ocos(j)b3 - 2.6 x IQ-'^^l^, 
+ 1.6 X 10~^ Ah'^Ar;.cos(f)b'r' + 1.1 X 10~^Ab'^ij,QCos(f)b'4 - 5.3 x lO'^A^^ 

- 2.6 X 10-'^At,Ab^cos(ptb - 2.7 x 10''^ At, AroCos(j)tr + 5.1 x 10~^AtoMiocos(j)n 
+ 2.3 X 10-%oM2oCOS(/)t2 + 4.1 x IQ-^Ao^socos^ta - 0.0016^1^^ 

+ 0. 00035 A^Miocos(l)ri + O.OOl^ro^2ocos0r2 - 3.7 x 10~^A^Msocos(l)r3 

- O.OOOSIA^, - O.OOO64A^^//oCos0^/4 + 0.00035^| 

+ 1.4 X 10-^ At'^Ab'^cos(l)t'b' + 1.4 X IQ-^^A.^cos^tv + 0.0002A^/^oCos(/)t/4 

+ 0.1571^1^0 - O.OOO13MioM2oCos0i2 + 6.4 x lO^^MioMaoCos^ia 

- O.OOllMlo + 2.9 X lO-^M2oM3ocos023 - 0.00019M|o 
+ 1.9 X 10-^m^^o - 0.0016m|^o - 1-3 x 10-^m|^o 

- 0.0016mf^o + + 1-8 x IQ-^m^^o 

- 1.3 X 10-^m^^o + 0-0025//(,^ (79) 

ml{tz) = 0. 00012 Ay^At,,cos(l)b't + 1.7 x 10~^Ab'^A^,cos(l)b'r + 0. 000069 Ab'^Miocos(l)b'i 

+ O.OOOSAb'^M2ocos4)b'2 + 0.0036^6^ M3oCos0fe/3 + 1.5 x 10"^^^^ A;,,, cos^^/^ 

- 1.2 X 10~'^ A^>AtoCos(j)r't - 0. 00052 A^i^A^^cos(()r'r + 0.000054^^/ Miocos(^^/i 
+ 0. 00015 A^>M2ocos(f)^,2 - 2.4 x 10~M^/M3ocos(/)^/3 - 0.00017 At'^AboCos(j)t'b 

- 0. 27 At>At„cos(t)t,t + 1.9 x 10-Mt^A,„cos</)tv + O.OlSAj/ Miocos0i/i 



+ Om2A^M2ocos(^t'2 + 0.39At'^M3ocos^t'3 - O.OSIM; 
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- 0.51M^Q + 0.96m^o - O.OOO44/i[,AboCos046 

- O.O98/ioA„cos04t - O.OOO24/ioA^„cos04r + O.OO79/XoMiocos04i 

+ O.O52/i[,M2oCos042 + O.26/i[,M3oCos043, (80) 

B.2 High tan/3 regime 

m]j^{tz) = O.OlAAl^ -0.0012Ab,AroCos(t)br- 0.0017 Ab^Mwcos(t)bi 



2 







- 0.014^6o^20COS(/)62 - O.O65A6oM3oCOS063 - 0.18^^, 

+ 0.0A4.AbiA^^cos(t)bir' - 0. 035 Ay^^cosctib' A - 0.083A^^ 

+ O.OlAt^Ab^cos^th - O.OOO53Ao^roCos0t^ + O.OlAoMioCos^a 

+ O.O6Ao^2ocos0t2 + O.27AtoM3ocos0t3 - 0.0011< 

+ O.OOO28A^oMioCos0^i + O.OO14A^oM2oCos0^2 + O.OO49A^oM3ocos0^3 

- 0.056v4^, - 0.03lA^^//oCOS(/)^/,4 + 0.096^| 

- 0. 032 At'^Ay^cos^)^^ + 0.0049A;, ^r^cos^t/,-/ + O.O3Aj,A*oCos0t'4 
+ 0.001371^1^0 - O.OO5MioM2ocos0i2 - 0.025MioM3oCOS(/)i3 

+ O.2M20 - O.17M2oM3ocos023 - 2.6M|o 
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+ 0.029m^^o + 0.028m^^o + 0.6m^„o 

- O.OOiemf^o - 0.0016mf^o - 0.37ml 

- O.W - 0.0083/x(f , 



= -O.llAl^ + 0m3Ab,Ar^cos(l)bT + 0.0025A„Miocos(/)w 

+ 0.069Ab^M2ocos(f)b2 + O.SGAboMaocos^ba + O.OSIA^^ 

- O.OO74A5/ A^/cos06v/ - O.OOO43A6^//[,cos0y4 + 0.009^^^^ 

+ 0.02lAt,Ab,cos(f)tb - 0.0032 At, A^^cos(l)tr - 0.0013A„Miocos(/)a 

- O.OUAt,M2oCos(l)t2 - 0.069AoM3oCOS(/)t3 - 0.046^2^ 

+ 0.0096AroMioCOS(l)rl + 0.018AroM2oCOS(l)r2 - 0.053AroM3oCOS(l)r3 

+ 0.01 lA^, + O.OO45>l^///[)Cos0^/4 - 0.24^^^^ 

- 0. 025 At'^Ab'^cos(f)t'b' + 0. 001 At'^ATi^cos(l)t'r' - 0.1lAt'^HQCos(f)t'4 
+ O.OllMi^o - O.OO5MioM2ocos0i2 - 0.0055MioM3ocos(/)i3 

+ 0.22M|o - O.I6M2OM3OCOS023 - 2.1M3^o 

- 0.31m^^o + 0.61m^^o + 0.03m^^o 

- 0.077mf^o - 0.077mf^o - 0.28ml, 
+ 0.03m,2 . + 0.13/x^,^ 



= -0.036^1^0 + O.OOAAboA^^cos(f)br + 0. 0013 Ab^Mwcoscpbi 

+ 0.022AboM2ocos(f)b2 + 0.1 AboM3ocos(f)b3 - 0.041A^^ 

+ 0.0048^/^ A^/cos06V' - O.O15A^//oCos06/4 - 0.024^^^ 

+ 0.01lAtoAboCos(l)tb - O.OOO83At(,A^(,cos0t,- + O.OO28AtoMiocos0ti 

+ 0.015At,M2ocos(j)t2 + 0.067At,M3ocos(l)t3 + 0.0046A^^ 

- 0.00095AroMiocos(l)ri - O.OO42A^(,M2ocos0T-2 - O.OllA^oM3ocos0^3 
+ 0.0065^^, + O.OO46^^^//oCos0^/4 - 0.052^^^^^ 

- O.OigAt'^Ab'^coscI)^^ + O.OOUAt'^Ari^cos(l)t'r' - 0.029At'^ij,'oCos(l)t'4 

- O.OllMi^o - O.OO19MioM2ocos0i2 - O.OllMioM3ocos0i3 
+ 0.32M22o - O.IIM2OM3OCOS023 + 4.7M|o 

- 0.098m^^o-0.091m^^o-0.12m^^o 
+ 0.0072mf^o + 0.0072m;^^o + O.lSm^Q 

- 0.12ml, + 0.35i,'l 



= 0.0094^^^ - 0.00082 Abo A^ocos(f)br - 0.001 lA^Miocos^bi 

- 0.0095AboM2ocos(l)b2 - O.O43A(,M3ocos063 + 0.064^^^^ 

- 0. 01 Ab'^Ari^cos(l)b'r' + O.OO5lA^//oCos0b/4 - 0.055AI 
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+ 0. 0067 AtoAboCos(f)tb - 0. 00035 AtoAroCos(f)tr + O.OO66^t(,Miocos0ti 

+ OMAt,M2ocos(f)t2 + 0.18At^M3ocos(f)t3 - O.OOOTe^l^^ 

+ 0.00019A^oMiocos(/)^i + O.OOO94A^(,M2ocos0^2 + 0.0033A^(,M3ocos(/)^3 

+ 0.017^^, + O.OO73A^///oCos0^/4 - 0.17A| 

- 0. 013 At'^Ai/^coscpt'iy + 0.0002AAt'^AT-i^cos(f)t'T' - O.O78At^//oCos0t'4 
+ 0.04371^1^0 - O.OO33MioM2ocos0i2 - O.O17MioM3ocos0i3 

- 0.19M|o - O.llM2oM3oCos023 + 4.6M|o 
+ 0.02m^^o + 0.018m|^o-0.27m|^o 

- O.OOllmf^o - O.OOllm^^o - 0.25mt\o 

+ 0.73m,2 . + 0.43/i(,2 (84) 



ml^itz) = -0.08lAl^ + O.OOS9AboAroCos(f)br + 0.0038AboMiocos(l)bi 

+ O.O53^6oM2oCos062 + 0.25AboMsoCos(f)b3 - 0.15^^^ 

+ 0. 02 Ab'^A^^cos(j)b'r' - 0. 036 Ab'yQCos(j)b' 4 + 0.0064A(g 

+ 0. 015 AtgAboCos(f)tb - 0. 0013 AtoA^cos(f)tr - O.OOllAo-^iocos^a 

- O.OlAt,M2ocos(l)t2 - OM7At,Msocos(l)t3 + O.OIA^^ 

- O.OO21^^oMiocos0^i - O.OO93A^oM2ocos0^2 - O.O25^^oM3ocos0^3 

- 0.0037 A^, + 0.0019A^i^ij,'oCos(f)r'4 + 0.066A^^ 

- 0.025At'^Ab'^cos(f)t'b' + 0. 0026 At'^Ari^cos(j)t'r' + 0. 02 At'^/iQCoscpt' 4 
+ O.OlMi^o - 0.00056MioM2ocos(/)i2 - 0.0055MioM3ocos(/)i3 

- 0.1471^2^0 - O.IIM2OM3OCOS023 + 4.9M|o 
+ 0.78m^^o-0.2ml,^o + 0.021ml,„o 

+ 0.015mf^o + 0.015mf^o-0.2m?,o 

+ 0.021m^^o + 0-28//0 , (85) 



mf^itz) = 0.007 + 0.02Ab,AroCos(l)br - 0.0032Ab,Miocos(l)bi 

- 0.011AboM2ocos(l)b2 - 0. 0082 AboM3ocos(l)b3 - 0.0049A^^ 
+ 0.023^6^ A^^cos^ft/^/ + 0. 01 Ab'yocos(j)b' 4 - 0.0005AJ 

- 0. 00072 AtoAboCos(f)tb - O.OO13Ao^ToCos0t^ + O.OOO27AoMoCos0ti 
+ 0.001 lAtoM2ocos0t2 + 0.0015At,M3ocos(f)t3 - 0.062^^^ 

+ 0.013AroMioCOS(l)rl + O.O32A^M2oCOS0^2 - O.O15A^oM3oCOS0^3 

- 0.064^^, - OMArl^lJ,'oCOS(f)r'4 + O-Ol^A}^ 

+ 0.00067 At'^Ab'^cos(l)t'b' + 0.0016Ai^ A^^cos^t/^' + O.OO65At^/XoCos0t/4 

+ 0.021Mi^o - O.OO42MioM2ocos0i2 + 0.0028MioM3ocos(/)i3 

+ 0.43M|o + 0.011M2oM3oCOS(/)23 + 0.043M|o 
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+ O.OlTm^^o - 0.084m^^o - 0.0011ml,^o 
+ 0.9ml,- 0.1ml, + O.Omml, 

- 0.0011m,2^o + 0.14//(?, 

mlitz) = O.OUAl + 0.039AboAr,cos(l)hr - 0.0063AboMiocos(j)bi 

- 0.023AboM2oCos(f)b2 - O.OlGAboMgocos^fts - 0.0097A^, 
+ 0.04:5Ab'^A^;cos(f)b'r' + O.O21Afe/^/ioCos06/4 - O.OOIA^^ 

- O.OOUAtoAboCos(f)tb - 0. 0025 AtoA^^coscptT + 0. 00053 AtoMiocos(f)n 
+ O.OO22At„M2ocos0t2 + O.OO29A„M3ocos0t3 - 0.12^2^ 

+ O.O25A^oMioCOS0^i + 0.0QAAr,,M20COS(/)r2 - 0.03A^^M3oCOS(pr3 

- O.UAl, - 0. 081 Ar'U'nCOS(f)r'4 + O.OSSA^, 

+ 0. 0013 At'^Ab'^coscpt'b' + 0. 0032 At'^Ar;cos(f)t>r' + 0. 013 At'^fj.QCOs4>t'4 

+ 0.12M2q - O.OO83MioM2ocos0i2 + O.OO56MioM3oCos0i3 

- O.llMlo + O.O23M2oM3ocos023 + 0.087^3^0 

+ 0.034m2^o - 0.17m|^o - 0.0022m|„o 

- 0.2ml, + 0.8ml, + 0.031ml, 

- 0.0022m'^t^, + 0.27i^'^, 

ml{tz) = 0. 0052 Ab'^AtgCos(l)b't + 0. 024 Ab'^AroCos(l)b'r + O.OO35Ab^MioCos06/i 

+ 0.062Ab'^M2ocos(f)b'2 + 0.31^6^ M3ocos06/3 + 0. 022 A,i^AboCos(l)r'b 

- O.OOIGA^^ Aocos0r't - 0.062A^/ A^ocos0^v + 0.0058^1^^ Miocos0^/i 
+ 0.0098A^^M2ocos(t)r'2 - O.O37A^^M3ocos0^/3 - 0. 0057 At>^AboCos(l)t'b 

- 0. 21 At'^At^cos(t)t't + O.OO19A^,^roCos0tv + O.O12Aj,Miocos0t/i 
+ 0. 078 At'^M2ocos(f)t' 2 + O.35A^,-^3ocos0t/3 - 0.036Mio 

- 0.36M|o + 0.68m3o - 0. 015 n',AboCos(l)4b 

- 0.04:2ii',AtoCos(f)u - 0.017 ii',Aj.^cos(j)4r + 0.0065ii',Miocos(j)4i 

+ O.O37//oM2oCOS042 + O.UfJ,',M3oCOS(j)43. 

B.3 trilinear terms in the NHSSM 

At the low values of tanP: 

At{tz) = -0.00063^60 + 0.22Aio + 3.6 x 10"^^^^ - 0.029Mio - O.23M20 - 1.9M; 

Ab{tz) = 0.99Ab^ - 0.13Ato - 0.00079A^o - 0.033Mio - O.48M20 - 3M30 

Ar{tz) = -0.0032Abo + 0.00029Ato + A^^ - O.I6M10 - O.53M20 + O.OO5M30 

At'{tz) = 0.00061>l6^ - 2.8 X 10"%^ + 0.49Atj^ + 0.46/Xo 

Ayitz) = 0.63^6^ - 0.00044A^/ + O.UAt>^ + 0.19^l', 

A'itz) = -O.OOI8A5/ + 0.49A^^ - 0.00026Aj, + 0.47/x(,. 
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When tan(3 is high: 

Atitz) = -0.05Afeo+0.21Ao + 0.0045A^o- 0.027Mio-0.21M2o-1.8M3o 

Ab{tz) = 0.38A6„- 0.072^0 -0.055A^o- 0.0092Mio-0.25M2o- 2. IM30 

A,{tz) = -0.26A, +0.027A„ + 0.62A,„-0.11Mio-0.32M2o + 0.44M3o 

Atitz) = 0.082^/ - 0.0065A^^ + 0.42 A;, + 0.18/i[, 

Ab'itz) = 0.54A. - 0.069A^^ + 0.12^;, + 0.083/io 

Ar'itz) = -0.29 Ah'^ + 0.6A^/ - 0.036^^/^ + 0.4/x[,. (90) 

Note that for the same values of tanf3 one-loop MSSM results can be obtained from the 
NHSSM solutions via the appropriate transformations (see text for details). 

C MSSM & NHSSM under universality assumption 

For the sake of simphcity and completeness, we also provide the solutions using (4), both in 
the MSSM and NHSSM; mass^ and trilinear terms are presented in the following subsections. 

C.l MSSM under universal terms 

With the help of (4) for low tan(5 MSSM results are 





= -0.087A^ + O.38A0M - 2.8M2 - 0.16mg, 




= -0.0033^2 + O.OllAoM + 0.49M2 + 0.99m^, 


mlitz) 


= -0.03^2 + 0.13AoM + 5.7M^ + 0.61m^, 


mlitz) 


= -0.058A^ + O.25A0M + 4.1M2 + 0.22m^, 


ml^itz) 


= -0.0017^2 + 0.0072AoM + 6.3M^ + 0.99m2, 


mlitz) 


= -0.00078^2 + 0.00067AoM + 0.52M^ + ml, 


mlitz) 


= -0.0016A^ + 0.0013AoM + 0.15M^ + m^, 


mlitz) 


= -O.SSAofiQ + OMmlo + 0.26M/io, 


Atitz) 


= 0.22^0 -2.2M, 


Abitz) 


= 0.074^0 - 0.3M, 


Aritz) 


= 0.052Ao - 0.036M, 



for high tan(3 MSSM results can be written as 

mjjjtz) = -0.061^2 + O.27A0M - - 0.12mg, 

mj^^itz) = -O.IA^ + 0.32^M - 2.M2 - 0.066m^, 

mlitz) = -0.04lA^ + 0.19^oM + 4.9M2 + 0.36m^, 

mlitz) = -0.041^^ + 0.18AoM + 4.3M^ + 0.25m^, 

mlJtz) = -0.042^^ + O.2IA0M + 4.7M^ + 0.46m^, 
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mf^itz) = -0.037^2 + O.OOQMoM + O.SIM^ + 0.75m^, 

mf^{tz) = -0.075^1^ + 0.02^0^ + 0.12M2 + 0.49mg, 

ml(tz) = -0.5Ao/Xo + O.eSmgo + 0.59M//0, 

A(i^) = 0.16^0 -2M, 

= 0.2lAo-2M, 

Ar(tz) = 0.2Ao + 0.0041M. (92) 



C.2 NHSSM under universal terms 

with the help of (4) again for low tan/S mass^ terms: 

mjf^{tz) = -0.087^2 + 0.10A[? - 0.16ml - 2MM^ + OMlA'^iio + 0.14^[f + O.SSAqM, 

mjj^ (tz) = -0.0033^2 - 0.37^0^ + 0.99m^ + 0.49 - 0.31 A^q + 0.6fi'^ + 0.011 AqM, 

ml^ (tz) = -0.03^2 - 0.089A[f + 0.61m^ + 5.7AP - O-OSA^ + O.Sfi'^ + O.ISAqM, 

ml^{tz) = -0.058^2 - 0.18 A'^ + 0.22ml + A.lK'f - 0.16A[,^o + 0.6/i(f + 0.2h A^M, 

ml^ [tz] = -0.0017^2 - 0.00079A[f + 0.99m^ + 6.3M^ - O.OOOISAX + 0.0029/^0^ 

+ 0.0072AoM, 

mf^ (tz) = -0.00078A^ - 0.00023A[f + ml + 0.52M^ - 0.00022yl>[, + 0.0012/x(f 

+ 0.00067^oM, 

mf^{tz) = -0.0016^^ - 0.00045A[f + ml + O.ISM^ - 0.00044AX + 0.0025//(f 

+ 0.0013AoM, 

ml{tz) = -0.27AoA() - 0.56M^ + 0.96m^o - 0-099AoAto + O.SAqM + 0.32/x[,M, 

At{tz) = 0.22Ao-2.2M, 

Ah{tz) = 0.86^0 -3.6M, 

Ar{tz) = 0.99^0 - 0.68M, 

At,(tz) = 0.49A(, + 0.46//(„ 

Ab'itz) = 0.77A(, + 0.19a('o, 

Ar'itz) = 0.49 + 0.47//;,. (93) 
For high tan(3: 

m]j^{tz) = -O.OeiAg - 0.12>l(f - 0.12mg - 2.6M2 - 0.036AX - 0.0083//(f + 0.27AoM, 

m]j^{tz) = -0.lA2-0.21A[f- 0.066m2-2.M2-0.1lAX + 0.13Ato +0.32AoM, 

(iz) = -0.041^2 - 0.1 A(f + 0.36m^ + 4.9M2 - O.OAA'^n'^ + 0.35/x[f + 0.19^oM, 

m2^(tz) = -0.041^2 - 0.1lA[f + 0.25m2 + 4.3M2 - 0.065AX + 0.43/i[f + O.ISAqM, 

{tz) = -0.042A^ - 0.086A[f + 0.46m2 + 4.7M^ - 0.014Ao/io + 0.28/x[f + 0.21 AqM, 

mj^ {tz) = -0.037A^ - 0.027A(f + 0.75mg + O.SIM^ - 0.023Ao//o + 0.14//(,^ 

+ 0.0099AoM, 
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mf^{tz) = -0.075^2 - O.OBAA'^ + 0.49m^ + O.llM^ - 0.047A>[, + 0.27^'^ + OmAoM, 

ml{tz) = -0.23^>1() - + 0.68m^o - 0.074Ao//o + O-SA'f^M + 0.19//[,M, 

= O.I6A0-2.IM, 

Ab{tz) = 0.25^0 -2.3M 

A^(tz) = O.39A0 + 0.0081M, 

Afitz) = 0.5>1(, + 0.18//(), 

Ab'itz) = 0.6A[, + 0.083/i[„ 

>lr'(iz) = 0.28A'o + 0.4/i[,. (94) 
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